Local unitary equivalence and entanglement of multipartite pure states 
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The necessary and sufficient conditions for the equivalence of arbitrary n-qubit pure quantum 
states under Local Unitary (LU) operations derived in [B. Kraus Phys. Rev. Lett. 104, 020504 
(2010)] are used to determine the different LU-equivalence classes of up to five-qubit states. Due 
to this classification new parameters characterizing multipartite entanglement are found and their 
physical interpretation is given. Moreover, the method is used to derive examples of two n-qubit 
states (with n > 2 arbitrary) which have the properties that all the entropies of any subsystem 
coincide, however, the states are neither LU-equivalent nor can be mapped into each other by 
general local operations and classical communication. 



I. INTRODUCTION 

The subtle properties of multipartite entangled states 
allow for many fascinating applications of quantum in- 
formation, like one-way quantum computing, quantum 
error correction, and quantum secret sharing [1, 2]. The 
theory of many-body states plays also an important role 
in other fields of physics which deal with many-body sys- 
tems [3]. Thus, the investigation of the non-local prop- 
erties of quantum states is at the heart of quantum infor- 
mation theory. Compared to the bipartite case, which is 
well understood, the multipartite case is much more com- 
plex due to the exponential grows in the dimension of the 
Hilbert space. Despite its relevance and the enormous ef- 
fort of theorists, many problems regarding multipartite 
entanglement are still unsolved [4] . Several entanglement 
measures for multipartite states, like the tangle [5], the 
Schmidt measure [6] , the localizable entanglement [7] , or 
geometric measure of entanglement [8] have been intro- 
duced. Moreover, different classes of entangled states 
have been identified [9] , and a normal form of multipar- 
tite states has been presented [10]. However, even for the 
simplest case of three qubits the entanglement proper- 
ties are still not completely understood. One of the main 
reasons for the lack of knowledge is, arguably, that for 
many-body entangled states we do have only few appli- 
cations [4] . This results into the existence of few known 
operational entanglement measures. 

One approach to gain insight into the complicated 
structure of multipartite states is to consider a re- 
stricted class of states, like for instance stabilizer states 
[1], matrix-product states [11], projected entangled 
pair states [12], Locally Maximally Entangleable States 
(LMESs) [13], or Gaussian state [14]. Considering a re- 
stricted set of states enabled researchers to gain a lot of 
intuition about the usefulness and manipulation of them. 
This knowledge in turn led to many of the fascinating ap- 
plications of multipartite states. 

Another way to gain insight into the entanglement 
properties of quantum states is to consider their intercon- 
vertability. That is, given two states \^f), |<I>) the question 
is whether or not I*) can be transformed into |<I>) by lo- 
cal operations [4] . One particularly interesting case is the 
LU-equivalence of multipartite states. We say that a n- 



partite state, |\P) is LU-equivalent to |4>) (|^) ~ LU |$)) 
if there exist local unitary operators, U\,...,U n , such 
that [*) = U\ ■ ■ ■ ® U n |$). Note that two states which 
are LU-equivalent are equally useful for any kind of ap- 
plication and they posses precisely the same amount of 
entanglement. Another insight is gained by considering 
more general operations, like (deterministic) local oper- 
ations and classical communication (LOCC). Since the 
implementation of such operations does not consume en- 
tanglement, a state, |^) which can be mapped into the 
state |$) by LOCC is necessarily at least as entangled as 
|$). Thus, all these investigations of convertibility lead 
to a new insight into general problem of characterizing 
the different types of entangled quantum states. 

In this article we will mainly focus on the LU- 
equivalence of multipartite states. Local polynomial in- 
variants have been introduced to distinguish the differ- 
ent LU-equivalence classes [15]. However, even though 
it is known that it is sufficient to consider only a finite 
set of them, this complete finite set is known only for 
very few simple cases. In [16] a method to solve the 
LU-equivalence problem for arbitrary n-qubit states has 
been presented. There, an algorithm which determines 
the local unitaries, which map the states into each other 
(if they exist) has been derived. Within this algorithm, 
different classes of states, which are easily characterized, 
are distinguished. It has been shown that two states 
which arc within two different classes cannot have the 
same entanglement. 

Here, we will use the algorithm presented in [16] in 
order to investigate the non-local properties of multipar- 
tite states. We will present the LU-equivalence classes of 
few-body systems and obtain a new insight into multi- 
partite entanglement. The main results derived here will 
be summarized in Sec II. 

The sequel of the paper is organized as follows. After 
presenting the main results of this article (Sec II), we 
review the necessary and sufficient conditions for LU- 
equivalence derived in [16]. In Section IV we will derive 
some additional methods to determine the local unitaries 
(if they exist) which interconvert the two states. In Sec- 
tion V we will characterize the LU-equivalence classes of 
up to 4 qubits. For five-qubit states we consider the most 
challenging class (for using the algorithm) and show how 



the local unitaries can be determined then. For an arbi- 
trary n (with n > 2) the existence of n-qubit states, \^f) 
which are not LU-equivalent to their complex conjugate 
will be shown by presenting examples in Section VI. In 
Section VII it will be shown how the algorithm can be 
employed to solve the LU-equivalence problem for cer- 
tain mixed states and also states which describe d-level 
systems. The new insight gained into multipartite entan- 
glement will be discussed in Section VIII. 

Throughout this paper the following notation is used. 
The Pauli operators will mainly be denoted by X,Y,Z. 
Whenever we need the whole set of Pauli operators we 
will use the notation Si = X, S2 = Y, and , S3 = Z 
and H denotes the Hadamard transformation. Other- 
wise, the subscript of an operator will always denote 
the system it is acting on, or the system it is describ- 
ing. The reduced states of system i\,...ik of |\fr) (\&)) 
will always be denoted by Pi 1 ...i k (o'i 1 ...i k ) resp., i.e. 



Pii...ik 



tr- 



, (|*) ($|). We denote by i the classical 



,i n ) with i k <s {0,1} Vfc e {1, ...,n} 
.,i n ) denotes the computational basis. 



bit-string (ii, . 
and |i) = |ii, . 
Normalization factors as well as the tensor product sym- 
bol will be omitted whenever it does not cause any confu- 
sion and 1 will denote the normalized identity operator. 
The eigenvalues of some matrix M will be denoted by 
eig(M). For a subsystem A we will denote by Ea(\*&)) 
the bipartite entanglement between A and the remain- 
ing systems measured with the Von Neumann entropy 
of the reduced state, pa- For instance, Ei(\^>)) = S(pi) 
will denote the entanglement between qubit i and the 
remaining n — 1 qubits. As commonly used, the states 
|$±) = |00) ± |11), |* ± ) = |01) ± |10) denote the Bell 
basis. The state |\&*) will denote the complex conjugate 
of the state |\I/) in the computational basis and M T will 
denote the transpose of an operator M in the computa- 
tional basis. 



II. MAIN RESULTS 

On the one hand, the criterion for LU-equivalence [16] 
will be used to characterize the LU-equivalence classes 
for few-body states. On the other hand, it will be em- 
ployed to shed new light on multipartite entanglement. 
Furthermore, the criterion will be generalized to certain 
mixed states and also states which describe d— level sys- 
tems. The main results derived here can be summarized 
as follows. 

i) Characterization of LU-equivalence classes: The 
LU-equivalence classes of quantum states describing up 
to five qubits will be characterized. For two, three and 
four qubit states all the classes will be explicitly derived. 
For five qubit states, whose classification would work 
analogously, only the most challenging subset of states 
will be considered. It will be explicitly shown how the 
algorithm can be used to determine the local unitaries (if 
they exist) which transform one state into the other. 

ii) New insight into multipartite entanglement: 



a) The algorithm presented in [16] distinguishes be- 
tween different classes of states, like the class of states 
with pi2 = pi (8> I2 and the one with pi2 ^ pi ® 12- It 
will be shown how this classification enables us to gain a 
new insight into multipartite entanglement. For instance, 
one class would be the one where one of the two-qubit 
reduced states is completely mixed. For four-qubit states 
it will be shown that this class is completely character- 
ized by only three non-local parameters. Moreover, it 
will be proven that two states within this class are LU- 
equivalent iff the corresponding sets of three parameters 
coincide. Naturally, all the entanglement contained in a 
state within this class is also determined by those three 
parameters to which the following operational meaning 
can be given. Recall that the completely positive map 
(CPM), £<&, corresponding to a state |W) via the Choi- 
Jamiolkowski isomorphism [17] can be implemented us- 
ing a system prepared in the state |\I/) and local oper- 
ations. It will be shown that the non-local content of 
the CPM, £<$, is characterized by the three parameters 
mentioned above and vice versa. This leads to the new 
approach of characterizing the entanglement of a multi- 
partite state by the entangling capability of the operation 
which can be implemented using the state as the only 
non-local resource. 

This suggests a new method to characterize the en- 
tanglement contained in an arbitrary multipartite state: 
First divide the Hilbert space into the entanglement 
classes resulting from the algorithm in [16]. Note that 
these classes can be easily characterized. Then, the en- 
tanglement of a state within a certain class should be 
qualified and quantified. Probably, the different classes 
might also lead to different applications. For instance, for 
error correction, one way quantum computing and quan- 
tum secret sharing, we have that all the employed states 
have the property that all single qubit reduced states are 
completely mixed. 

b) The other new insight into multipartite entangle- 
ment which we will derive here using the LU-equivalence 
criterion is the following. For any n > 2 examples 
of n-qubit states, \^>) and |$) which have the proper- 
ties that for any subsystem A, composed out of arbi- 
trary many qubits, the eigenvalues of the reduced states, 
p A = tr^(|*) (*|) and a A = tr^(|$) ($|) coincide, 
will be presented. Therefore, all the bipartite entangle- 
ment in those two states, measured with the von Neu- 
mann entropy of the reduced states, coincide. It will 
be shown, however, that the states are neither LU- 
cquivalcnt nor LOCC comparable. Therefore, neither |\&) 
can be mapped into |$) by LOCC nor vice versa. Surpris- 
ingly, in those examples we will have |$) = |^*), where 
\^>*) denotes the complex conjugate of |\I/) in the com- 
putational basis. The fact that \§>) and |\&*) can have so 
different non-local properties does not seem very phys- 
ical. As a consequence of the existence of these states 
it will be suggested to divide the Hilbert space into two 
subsets, in case |\1/) is not LU-equivalent to [$*). One 
which corresponds to |\&) and one which corresponds to 



its complex conjugate. The non-local properties should 
then be investigated within one of the subsets since there 
will not be a physical measure which will distinguish be- 
tween |*) and |**). 

iii) Generalization of LU-equivalence criterion: It will 
be demonstrated how the solution of the LU-equivalence 
for pure n-qubit states can be generalized to mixed states 
and also to d-level systems. 



III. LU-EQUIVALENCE OF MULTIPARTITE 

STATES 



Here, we briefly summarize the necessary and suffi- 
cient conditions for the existence of LU operations which 
transform two n-qubit states into each other [16]. We 
will first review a standard form for multipartite states 
(see also [18] and [13]) and will provide some examples 
of states in their standard form. It has been shown that 
two generic multipartite states, i.e. states where none 
of the single qubit reduced states is proportional to the 
identity, are, like in the bipartite case, LU-equivalent iff 
their standard forms coincide [16]. For non-generic states 
the systematic method to determine the local unitarics 
(if they exist) which interconvert two arbitrary states will 
be reviewed. 



A. Standard form of multipartite states 

Let us first recall the definition of the standard form for 
multipartite states. Any decomposition of a multipartite 
state which has the property that the single qubit re- 
duced states are all diagonal in the computational basis 
is called trace decomposition. It is obtained by applying 
local unitary transformations, £//, which diagonalize the 
single qubit reduced states, pi, i.e. U} Pi{Ul)^ — Di = 
diag(A|, Af ). A sorted trace decomposition, which we will 
denote by \^ s t) in the following, is then defined as a trace 
decomposition with A* > Af . The sorted trace decompo- 
sition of a generic state, | "J) with pi ^ \ Vz is unique up to 
local phase gates. That is U\ . . . U n \^ s t) is a sorted trace 
decomposition of a generic state, | "J/) , iff (up to a global 
phase, ao) Ui — Zi(ai) = diag(l,e Mi ). It is straight- 
forward to impose certain conditions on the phases «», 
i € {0, . . . , n} in order to make the sorted trace decom- 
position of generic states unique [16]. We call this unique 
sorted trace decomposition standard form of the multi- 
partite state. Note that any state can be transformed 
into its standard form by local unitary operations. 

Let us now also recall how the standard form can 
be defined for states with pi = 1, for some system i 
[13]. In this case the standard form can be chosen to 
be lim e ^ |*( e )), where | \I/ (e)} denotes the unique stan- 
dard form of \/\ — e \^) + \/e|0), where the phase gates 
are fixed by the same conditions as for generic states [13]. 
It should be noted here that for non-generic states the 
standard form is not unique, as can be seen by the fol- 



lowing simple example of a three-qubit states. Both, the 
GHZ-state, |*) = |000) + |111) and HHH\^>) are stan- 
dard forms of the state |\1/), however they do not coincide. 
Let us now derive the standard form for some exam- 
ples. For 2-qubit states the standard form coincides with 
the Schmidt decomposition [19]. In [20] the standard 
form of three-qubit states has been derived. In order to 
present the standard forms of certain n-partite states we 
recall here the notion of the so-called Locally Maximally 
Entangleable States (LMESs) [13]. LMESs have been in- 
troduced as a new, physically motivated, classification 
of pure quantum states describing n qubits. A state is 
called LMES if local auxiliary qubits can be attached to 
the system qubits in such a way that the resulting state 
is maximally entangled in the bipartite splitting system 
qubits versus auxiliary qubits. To be more precise, a 
state |\&) is a LMES if there exist local control operations 
Ci = |0) (0| U? + |1) (1| U}, with C/° 4 single qubit 
unitary operators acting on the system qubit i, such that 
the 2n-qubit state d C 2 . . . <g> C n |*> |+)® n , with 
|+) = l/\/2(|0) + |1)), is a maximally entangled state be- 
tween the system and the auxiliary systems. This set of 
states coincides with the set of states which can be used 
to encode locally the maximum amount of n independent 
bits. Prominent examples of these states are the stabi- 
lizer states, which are used for quantum error correction 
and one-way quantum computing. In [13] it has been 
shown that a state is LME iff it is LU-equivalent to a 
state of the form 



I*) 



>n / t 



e -«|i) = J7*|+) 



(1) 



where cti € R and UK denotes the diagonal unitary op- 
erator with entries e Mi . Thus, a state is LME iff there 
exists a product basis such that all the coefficients of the 
state in this basis are phases. 

Note that all those states can be easily transformed 
into their trace decomposition by applying the local uni- 
tary operations HUi, where Ui = diag(e I( ^, 1), with 

cot(<fc) = 4f4 if (Yi) ^ and & = else. To de- 
rive from the trace decomposition the standard form one 
simply has to impose the conditions on the local phase 
gates, as mentioned above. 



B. Criterion for LU— equivalence 

Since the standard form is unique for generic states 
we have, similarly to the bipartite case that two generic 
states are LU-equivalent iff their standard forms are 
equivalent. 

Let us now turn to the more complicated case of non- 
generic states. First, the condition of LU-equivalence for 
generic states is rewritten in the following way. It can be 
easily seen that the standard forms of two generic states, 
|\&) and |$) are equivalent iff there exists a bit string 



k = (fci, . . . , k n ), local phase gates Zj(a,), and a global 
phase a s.t. 



•Z i {a i )X* i W i \%) = (>Z\V i \$) 



(2) 



where Wi (Vi) are local unitaries which diagonalize pi 
(ai). That is (g^ Wi |W) and ® 4 Vi |3>) are trace decom- 
positions of |^) and |$) resp.. For generic states, fcj is 
chosen such that the order of the eigenvalues of the sin- 
gle qubit reduced states of (g). X^Wi |*) and ® 4 Vi |$) 
coincides and the phases on are chosen to fulfill the con- 
ditions mentioned above [16]. Note that the reason for 
the freedom of the phase gates in Eq. (2) is simply due to 
the fact that we have been considering only single qubit 
reduced states to define the trace decomposition of mul- 
tipartite states. 

Obviously, two arbitrary states |\t) and |<&) are LU- 
equivalent iff there exist local unitaries Vi and Wi and a 
bit string k = (k±, . . . , fc„), and phases on s.t. Eq (2) is 
fulfilled. For non-generic states, a constructive method 
to determine the unitaries Vi and Wi in Eq. (2) has been 
presented in [16]. Once those unitaries are fixed it is then 
easy to decide whether or not there exist local phase gates 
for a certain bit string k such that Eq. (2) is fulfilled (see 
Lemma 1 below). 

Since we arc going to determine the local unitaries 
which transform two states into each other in Sec V 
we review here the constructive method to compute the 
unitaries V k ,W k . First of all, it is easy to see that if 
\ty) is such that there exists some system i such that 
Pi 7^ 11 the unitaries Vi and Wi can be determined by 
considering the necessary condition for LU-cquivalence, 
Pi = Ui<JiU\. Analogously to the generic case, the equa- 
tion Di — diag(A 4 1 ,A|) = WiPiWJ — Vi<JiV} determines 
Wi and Vi (and fcj = 0) uniquely up to a phase gate. 
Thus, for this case we have that |vE') ~ LU |$) iff there 
exist two phases, on and ao and local unitaries Uj such 
that 

i (l\ W^ s ) = e i(aa+a ^ (g) U 3 i (l\ Vi$ s ) for I E {0, 1},(3) 

where Wi, and Vi are chosen such that Di — 
diag(A^, \\) = WiPiW} — ViViV} . Hence, if there is one 
system where the reduced state is not completely mixed, 
then the problem of LU-equivalence of n-qubit states 
can be reduced to the one of (n — l)-qubit states. This 
statement can be easily generalized to the case where 
more than one single qubit reduced state is not com- 
pletely mixed. 

Let us now turn to the remaining case where pi = 
1 Vi. Instead of considering the necessary conditions 
Pi = UiCTiUJ , one considers the necessary conditions 

Pm,...,ni,k = U rll . . . Un l Uk<T ni ,...,n l ,kUl 1 ■ ■ ■ U^U^ for 

some appropriately chosen set {ni, . . .ni,k} and com- 
putes Uk as a function of U ni , . . . U ni . More precisely, 
it has been shown that if |\&) = U\ . . .U n \&) and if 
there exist systems m, ■ ■ .ni and k such that p ni .... ni .k 7^ 



Pm,...ni ® Ife then Vk in Eq. (2) can be determined from 
the state |$) and W k can be determined as a function of 
the unitaries U ni , . . . U ni . To be more specific, we assume 
without loss of generality that n\ = 1, . . . n; = I. Due 
to the condition pi t ...,i tk =h P1....1 Ife it can be shown 
that there exist at least two tuples i = (i\, . . . , i{) and 
j = (ji,...ji) such that at least one of the hermitian 
2x2 matrices B\ = A\ + (4J)t and C{ = iA\ - i{A\)\ 
where A\ = tr-,fc[|i) (j| |$) ($|], is not proportional to 
the identity. W. 1. o. g. we assume that 1 9^ B\ = 
tr^ fc [(|i) (j| +h.c) |$) (*|]. Using that \^>)=U 1 ... U n |$) 
we have 

U k B\ul = tr^ fc [(|i) (j| + h.c) ■U\...U}\*){*\Ui... Ui].{A) 

Since B\ is hermitian we can diagonalize it as well as 
the right hand side of Eq (4). It can then be shown 
that [*) = Ui . . . U„ |$) iff there exists i k <G {0, 1}, and 
a and a k such that e ia °X lk U(a k )W k (Ui, . . . ,Ui) |*> = 
U\ . . . V k ■ ■ ■ U n |$), where V k is the unitary which diago- 
nalizcs B? and can therefore be determined directly from 
the state |4>) and W k (Ui, . . . ,Ui) diagonalizes the right 
hand side of Eq. (4). 

Note that this constructive method to compute V k , W k 
is based on the necessary condition for LU-equivalence 
given in Eq. (4) for any ^-tuples i,j. Since the 2x2 
matrices occurring in this equation arc hermitian, one 
can, similarly to the previous cases, determine the uni- 
taries V k ,W k by diagonalizing these matrices. In con- 
trast to before we will find here, that W k might depend 
on Ui , . . . , Ui . Again, since those unitaries are obtained 
by diagonalizing a 2 x 2 matrix the phase gate occur- 
ring in Eq. (2) cannot be determined like that. This 
is the reason why the condition of LU-cquivalence has 
been rewritten in the seemingly more complicated form 
presented in Eq. (2). 

In order to check then whether or not there exist phases 
Oj such that Eq. (2) is satisfied, the following lemma 
[16], which will be proven here, has been used. We con- 
sider four n- qubit systems which will be denoted by 
A, B, C, D respectively. The i-th qubit of system A will 
be denoted by Ai, etc. Furthermore, we will use the 



notation \x), = (|0110) - [1001))^ 
E k Ik) (kk| 



and P l AC = 



Similarly, 



,Ci,Di 
,k I""/ \ w ^-\A 1 ,C 1 ,...Ai- 1 ,Ci- 1 ,A (+1 ,Ci +1 ...,A n ,C, 

we define P l BD for systems B,D. For a state \^) 
define Ky = {k such that (k| *) = 0} and 



m' 



and|* ) = |*) + 2X: 



-*Z^?=i' 



*{a,}) 



|k) for some phases c^ 



veK<s, 



Lemma 1. Let |\&) , |$) be n-qubit states. Then, there 
exist local phase gates, Zi(a.i) and a phase 01$ such 
that |*> = e ia ° (g)™ =1 Zi((Xi) |$) iff there exist phases 
{a*}"=o such that (i) I (i|*o)| = I (i|${a,})l Vi and (ii) 

(X^PXcPBD 



{1, 



'}• 



1*0^ l*o> B \^{ &1 }) C \^{ & ,}) D = Vz e 



Condition (ii) can be interpreted as follows. Taking 
two copies of the state |\&o) and two copies of the state 



^{Si}) an d projecting the four qubits, Ai, Bi, C\, Di onto where ag 



the state |x) i leads to a 4(n — l)-qubit state, which is in 
the kernel of Packed f° r an y system i. Before proving 
Lemma 1 we introduce here another lemma, which will 
be required for the proof. Using the same notation as 
before, we have 

Lemma 2. |\I/) can be converted into |$) by local unitary 
phase gates iff there exist phases {<5i}™ =0 such that \^o) 
is converted into I'&aJ by local unitary phase gates. 

Proof. If |*) = e iao 0™ =1 Z 4 (a,)|$), for some phases 
{at}, then Kq, = K$ and choosing on — ai for i £ 
{0, . . . , n} fulfills the condition. To prove the inverse 
direction we assume that there exist phases {<Xi}f =0 
such that l^o) = e ra ° <S)ILi ^i( a i) |^a«) f° r some phases 
{ai}. Due to the factor 2 in the definition of jv^o) 
and l^ai), this implies K^ — K$. Defining the pro- 
jector P = Ek^x* l k ) ( k l we have Pl^o) = |#) and 
Pe iao <S>7=i Z M)\^i) = e <ao ®2= 1 Z i (a i )P|# ai ) and 
therefore |*) = e ia ° (g)" =1 Z t (a t ) |$). 

D 

The reason for introducing this lemma is that it implies 
that if one wants to decide whether or not two states 
are up to local phase gates equivalent, one only needs 
to consider states where none of the coefficients in the 
computational basis vanish. Let us now use the lemma 
above to prove Lemma 1. 

Proof As mentioned above, due to the Lemma 2 
it remains to show that for any state \ip) with 
(k| V) + Vk we have that \ip) = e la « <g)? =1 Z t {a t ) \<f>), 
for some phases {a^ iff condition (i) and (ii) 
in Lemma 1 are satisfied. Note that condition 
(ii) is equivalent to (0k| V) (H| V>) <lk| <t>) <01| <j>) = 
(lk| tp) (01| \jj) (0k| 4>) (11| </>), where 0, 1 is acting on sys- 
tem i and k, 1 denote the computational basis states of 
the remaining n—\ qubits. 

Let us now prove the only if part: If \ip) = 
e la " 0" =1 Z(ai) \4>) for some phases {a,} then (i| ip) — 
e 1 ^'' {i\ (f>) , with <pi — ao + Efc afe * fe ' which implies 
(i). Condition (ii) (for i = 1) is then equiva- 
lent to e tWok+01l) x k i = e lWlk+0ol) x k i, where x k i = 
(0k| (p) (11| <f>) (lk| <p) (01| <p). It is easy to see that this 
condition is fulfilled since e l( * ok -* lk) = e~ Ml Vk. In the 
same way one can show that the conditions for i ^ 1 arc 
fulfilled. ' 

To prove the if part, we first note that condi- 
tion (i) implies that (i| v &) = e*^' (i| <fr), for some 
phases fa. Condition (ii) (for i = 1) implies 
then that e 1 ^^'^ 1 ^ = e^oi-^n) \/k, 1, since x ki = 
(0k| 0) (11| <j>) (lk| </>) (01| </>) ^ Vk, 1. Thus, e^"-* 1 ") 
must be independent of k and therefore, we have 



s J(0Ok-01k) 

have i 



e lai , for some phase a\. Equivalently, we 



J fe i 



e i(a 1 +0i k )^ w } iere qA > — _ ai anc [ q,( J _ 



.,(!) 



0. Similarly, we obtain e*^*! ** 

and therefore e"™!'**^---'*" = 
Continuing in this way we find e l ^ k 



k n — ( Pk 1 lk 3 



— „-4«2 



^(aP'+a^^+^n*., 



*°<g>r =1 ^ 



E< 



Thus, we have |-0) 



□ 



It is important to note here that the state on the right 
hand side of Eq. (2) is completely determined using the 
method summarized above. Thus, the set K^ in Lemma 
1 can be determined and therefore this lemma can be 
applied. The states are LU-equivalent iff the conditions 
in Lemma 1 are fulfilled for some bit string k. The uni- 
taries which interconvert the states are, up to the sym- 
metry of the states, uniquely determined and are given 
by Ui = W} Zi(ai)X ki Vi (up to a global phase) [32]. 

In summary, the LU-equivalcnce problem has been 
solved by presenting a systematic method to determine 
the local unitaries (if they exist) which interconvert the 
states. This has been achieved by determining Vi,Wi 
in Eq. (2) by imposing necessary conditions of LU- 
equivalence, like pi = UiOiU\ and Eq. (4). Once all 
the unitaries Vi, Wi are determined (even as functions of 
some others), the states are LU-equivalent iff there ex- 
ist local phase gates which interconvert the transformed 
states (after applying 0^1^, ® 4 Wi to |$), |*) respec- 
tively). This can then be easily decided by employing 
Lemma 1. 

Before ending this section let us present here another 
way of checking whether or not two states are intercon- 
vertible by local phase gates. Due to Lemma 2 we only 
need to consider states |\&), |3>) with Ky = K§ = 0. 
Here and in the following we will denote by Q the 
Hadamard product, i.e. the component-wise product 
and by /. we will denote the inverse operation, i.e. the 
component-wise division. For instance, if |^) = Ei a i I') 
and |$) = Ei & i|i>, with h ^ Vi, then |#)/.|$) = 
Ei 01/61 |i>- 

Lemma 3. Let |^) and |$) be n-qubit states with 
K^ = Kq, = 0. Then, there exist phases {ai} such 
that I*) = e ia ° (g^aii)!*) iff 1) |(i|*)| = |(i|$)| and 
ii) |\&) /. |<f>) is a product state. 

Proof. (Only if): If |#) = e ia ° ® 4 Z((Xi) |$) 
condition i) is obviously fulfilled. In order to 

show that condition ii) is fulfilled we use that 
e ia "<S>iZ{ai)\^) = e M0 (g) 4 Z(a J )|+>^ n OI$>. Thus, 



<a Z(on 



which is a product 



oiao Jllj a i k i 



l*)/.l#) = 

state. 

(If): Due to condition i) we have that |\I/) /. |<&) = 
5^e Mi |i), for some phases ai. That is, |\f) /. |$) is a 
LME state. Due to condition (ii) this LME state must 
be a product state, i.e ^e Mi |i) = (^ \<f>i), where 
\<f>i) = e J *o(A^ |0) + e l *iAi |1)) with X{ > 0. This implies 

that e iai = e^+^^y", where $ = Efe^o- Thus . 
the LME state is a product state iff it is equivalent to 
e ta ° <§Z) i Z{pn) |+)®" for some phases {af} and therefore 
|*) = e io "'® j ^(oi)|#). □ 

As mentioned above, Lemma 2 can be used to gen- 
eralize Lemma 3 to states, |W) for which Kq, 7^ 0. 
Note that condition ii) has a physical interpretation. 



The Hadamard product of two states \ip), \4>) corre- 
sponds to the state one would get by the following pro- 
cedure. Let \ip) {\<t>)) describe the system li,...,m 
(I2, ...,ri2) rcsp. and consider n pairs of maximally 
entangled two-qubit states, |<J> + ) = J2i=o N)> de- 
scribing systems I3, I4, . . .713,714. Then, \(j>}Q\ip) — 

®r=i«,i 2 ,i, iv') 1i ,...„ 1 i0)i 2 ,...„ 2 ®r=i i* + > i3li4 , where 

I "J ) denotes the GHZ states here. This resembles the 
procedure of gate teleportation [21]. Note that condition 
ii) is fulfilled iff there exists a product state, (^) i \<pi) such 
that|*o> = |*a 4 >O®il0i>. 



IV. ADDITIONAL METHODS TO COMPUTE 
THE LOCAL UNITARIES 

We have seen before how the local unitaries which oc- 
cur in Eq. (2) can be determined by imposing certain 
necessary conditions of LU-equi valence (see Eq. (4)). 
One might also use other necessary conditions for LU- 
equivalencc to determine those local unitaries. For in- 
stance, if I*) = Ui...U n \$) then tri(pi I2P12) = 

U 2 tri(<Jl l20"12)^2 and tr 23(/Ol23 ll'Pl'23 ll) = 

^i<g>6 r i/tr 2 3(cri23<8>li'cri/23<8>li)?/ ] ^>U{, . Of course, any 
generalization of these equations must be fulfilled too. 
Here we will use those and other necessary conditions 
for LU-equivalence to derive some additional methods 
to compute the unitaries in Eq. (2) for certain multipar- 
tite states. Depending on the properties of the states of 
interest one method or the other might be better suited. 
In Sec V we will use the various methods to compute 
the local unitaries directly, i.e. not as a function of 
other unitaries. This makes the characterization of LU — 
equivalence classes easier. 

Here, we will first consider the LU-equivalence of two- 
qubit mixed states. Then we will focus on those states 
for which there exists at least one system i with pi 7^ 11 
and will derive a simple way to determine the unitaries 
in Eq. (2). 



A. Two-Qubit Mixed States 

For two-qubit mixed states, p, o~, necessary and suffi- 
cient conditions for LU-equivalence have been derived in 
[22]. However, if p = pij (a = crij) denotes the reduced 
state of some systems i, j of a multipartite state, |\&) (|$)) 
resp. and the aim is to investigate the LU-equivalence of 
|\&) and |$), then one must determine all local unitaries, 



Ui,Uj, which fulfill pij 



UiUjCTijUjUj and then check 



if there exists one of them, which transforms the multi- 
partite states into each other. We are going to show here 
how to achieve this task. 

We have seen above that if there exists some sys- 
tem i such that pi 7^ 1, then Vi, Wi and fc, in Eq. 
(2) can be determined by imposing the necessary con- 
dition pi = UidiU} . Thus, it remains to consider the 



case where both reduced states are proportional to the 
identity which implies that p = 1 + Y] k , AjyEfc £;, 
where A = J2ki ^ki\k) (l\ is real. Applying the lo- 
cal unitary operation U\ U 2 to the state, p leads 
to Ui U 2 pU\ ®U\ = 1 + Efc,z A fc,i Sfe ^h with 
A' = OiAOj. Here, Oi,02 are real orthogonal matrices 
which are defined via the equation Ui(na)U] = (Oifta) 
for i = 1,2. Using the singular value decomposition of 
the real matrix A, A = 0\D0 2 ', where D is a diagonal 
matrix, and 0\^ are real and orthogonal, and the fact 
that the state 1 + Efc l Dk,k^k £fc, is Bell-diagonal 
shows that the eigenbasis of any two-qubit density ma- 
trix with completely mixed reduced states is maximally 
entangled. 

In order to show now under which conditions two two- 
qubit states are LU-equivalent we recall the following 
Lemma which was proven in [23] . 

Lemma 4. Any two maximally entangled basis of two 
qubits can be mapped into each other using local uni- 
tary operations [23]. That is, if {|*j)}f =1 and {|$»)}i=i 
denote two maximally entangled bases then there ex- 
ist four phases ji, and local unitaries U\,U 2 such that 
|#ij = e** J7i U 2 |$i) V« e {1, 2, 3, 4}. 

This lemma together with the fact that the eigenbasis 
of any two-qubit density matrix with completely mixed 
reduced states is maximally entangled implies the follow- 
ing corollary. 

Corollary 5. Let p, a be two-qubit density matrices 
with completely mixed reduced states. Then p ~lu o~ 
iff eig(p) = eig(a). 

Let us now consider two LU-equivalent states, p, a 
with pi = Cj = 1 for i = 1,2 and derive some con- 
ditions on the local unitary operations, which trans- 
form a into p. First we apply local unitaries, W{,Vi 
such that p = W{W 2 pWlwl = 1 + £;(£>,,);£;£; and 

a = V{V 2 aVM = 1 + £i(Ar)i£iEi, with D p = D a = 
diag(Ai, A2, A3). We choose w.l.o.g. the order of Aj such 
that if there is no degeneracy Ai > A2 > A3, else Ai = A 2 . 
If Dp is not proportional to the identity it is easy to 
see that p = U\U 2 aU{Ul implies that Ui is of the form 
Z(ai)X ki , for some phase Qj and fc» <G {0, 1} [33]. Thus, 
if p and a denote for instance the reduced state of system 
1 and 2 of some multipartite state, |\&), |$) respectively, 
then I*) ~ LU |$) iff Eq. (2) is fulfilled for V 1 ,V 2 and 
W ll W 2 such that W x W 2 pW\wl =0 = V x V 2 aV^ = 
1 + J2 t ( D p)i^i^ where D p = diag(Ai, A 2 , A3) with \ t 
sorted as mentioned above. 

Otherwise, if D p is proportional to the identity and p 7^ 
1, we apply the local unitaries Vi and Wi defined above 
and denote the resulting states again by p, a respectively. 
In this case we find p = a = 1— A l^ - ) {9~\ for some A 7^ 
0. Then, any pair of unitaries U\, U 2 which transforms 
a into p must fulfill that U 2 — U\. Hence, in this case 
we have that if p and a denote for instance the reduced 
state of system 1 and 2 of some multipartite state, |^), 



|$) respectively, then |\p) ~ LU |$) iff Eq. (2) is fulfilled 
for V\ = V 2 = 1, k x = k 2 = 0, Qi = a 2 , and Wi ~ W 2 — 
e J/3lXi e 47lZi , for some phases /3i,7i. Note that if pi 2 = 
1-A|* _ ) (*"| with A ^ 1 then |*) ~ LU |$) implies that 
i 2 (*~ |*) ~l;7 12 (*~ |$) since 17 J7 |*~) = |* _ > for 
any unitary U . Thus, similar to the case where p, 7^ 1 one 
would simply measure those systems where the reduced 
state is a full rank Werner state [31]. 

In the remaining case, that is if p = 1, it is clear that 
considering the two-qubit reduced state will not help us 
to find any condition on the local unitaries. 

So far we have seen that whenever there exist two 
systems i, j such that pij is not LU-equivalent to 1 + 
A J2i Si Si, then the unitaries Vi,Vj,Wi, Wj and ki, kj 
in Eq. (2) can be easily determined. We are going to 
show next that in this case also other unitaries, Vi and 
Wi, for I ^ {i,j} can be easily computed. As before 
we consider the case where pi = 1 Vi. In the following 
lemma we will say that the unitaries can be determined 
by considering a certain operator, if they can be deter- 
mined using the fact that the operator for the state \$) 
and the one for the state |$) must be LU-equivalent if 

Lemma 6. If there exist systems i,j such that pij is not 
LU-equivalent to 1 + X J^. Ej Ej, for any X £ R, then, 
for any system I for which either pu ^ 1 or pji =/= 1, 
Vi , Wi and ki can be determined by either considering pu 
or pji or by considering tri(pij pu) ortTj(pijPji). 

Proof. If pu or pji is not LU-equivalent to 1 + A J2i ^» ® 
Si, for some A, then V/,W and fc; can be determined 
as shown above. Otherwise, we assume without loss 
of generality that pu =/= 1. Then we have that pu = 
1 + ^2 t i _ 1 Ai 1 i 2 T,i 1 Ej 2 with A proportional to a 

real orthogonal matrix and p^ = 1 + ^ ; A^^E^ E; 2 
with A not proportional to a real orthogonal matrix. 



Then, we find tr i(pijPu) 



Z2i lt i 2 Aj 1 Z 2 Aj 1 j 2 E; 2 Ej 2 , 



where the matrix (A) T A is not orthogonal. Since 
the unitaries Wj,Vj are already fixed, the equation 

WjWitn{ P ijpu)W}w} = VjVfriipijOifiVJV? deter- 
mines Vi,Wi and fc; = 0. □ 



B. States where there exists a system i with pi ¥ 1 

Let us now turn to the case where there exists at least 
one system i such that its reduced state is not completely 
mixed. Without loss of generality we chose i = 1. Let 
us assume that the states of interest, |^) and |$) do 
have sorted trace decomposition. That is, in particular 
the unitaries V\ , W\ which make p\ and o~\ diagonal in 
the computational basis (see Eq. (2)) have already been 
applied. Then we know that U\ = Z(a\) for some phase 
a\. We will present now several methods to determine 
the unitaries Vj, Wj of Eq (2). As in [16] the idea is to 
construct out of the states |\I/) and |$) non-degenerate 
2x2 matrices which must be LU-equivalent in case \ty) 



and |$) are, e.g. UpiW = Oj, or tri[(pi l»)pi»] = 
L/jtri \{<j\ \)aii]U}. Those necessary conditions of LU- 
equivalence can then be used to fix Wi,Vi and fcj in Eq. 
(2). The local phase gates, which cannot be fixed in this 
way, must be determined at the end using one of the 
lemmata in Sec. III. 

Since the states \*f?) and |<J>) have sorted trace decom- 
position, we have |\&) = Jp\ 
with (^il^,) = Si.j and p\ 



^|0)|$ ) + \/T^|l)l$o> 



l*o) + Vi 


-Pi 


|1>I*1>, 


> 1/2 


and 


|$) = 


vith (<f>i $. 


/> = 


6ij and 



p 2 > 1/2, which is just the Schmidt decomposition for 
the bipartite splitting, system 1 versus the rest. Then we 
have that the states are LU-equivalent iff 1) p\ = p 2 and 
2) there exist phases, <f>, ot\ and unitaries Uj such that 



(5) 



j¥1 



where 71 = <f> + a.\ and 72 = 4> — ct\ (see Eq. 
(3)). Note that the last two conditions are fulfilled 
iff /3i |*o) (*o| + h |*i) (*i| - ® j¥1 UjifSx |$ > (*o| + 
P2 |*i) (*i|) ® j¥1 U] for all values of ft, ft. 

There are several ways now to compute the unitaries 
Vi, Wi in Eq. (2). First of all, if pi ^ 1 then U t = Z{a t ) 
for some phase on. Let us now consider the case where 
Pi = 1. If pu / pi ® 1 Eq. (4) can be used to compute 
Vi,Wi, and fcj. In particular, we would consider one of 
the matrices B^ m = tr^(|*j) (* m | + h.c.) or C Lm = 
tr^i(i \i&i) (ty m \ + ft.c.) for some I, m and diagonalize this 
matrix in order to compute Wj, Vj and ki- 
lt is the aim of this section to derive some other meth- 
ods to determine those unitaries. First, we will use the 
fact that Eq. (4) must be fulfilled for any values of I, m if 
the states are LU-equivalent. Considering certain combi- 
nations of those equations will lead to other approaches 
to determine the unitaries Vj, Wi and the bit value ki. In 
the second part of this section we will show how a com- 
bination of the Eqs. (5) can be used to compute those 
unitaries. 

If pij 7^ p\ 1 then Vi and Wj can be easily computed 
as follows. First of all, it is clear that if p\ 7^ u\ then 
they states are not LU-equivalent. Thus, we assume that 
p\ = o~\. The fact that p\ 7^ l and pi = \ implies that 

Pii = l + aZl + J2l d2=l ^ jl 
and where A / since pij 7^ p\ 
<7ij = 1 + aZl 



1 J2 Sji 



T,j 2 , for some a € R 



E 3 



1. Similarly we have 

,., , ;j =i r ji h S ii ® S J2 • As mentioned 
before, the two states are LU-equivalent, i.e. pij = 
U\UiO-\iU\Ul iff there exists a real orthogonal matrix Oj 
and a phase oc\ such that A = OiTOf , where 0\ — 
R z (a\). As explained in Sec. II A, if A is not orthog- 
onal, then the unitaries Vj, Wi can be easily computed. 
Otherwise, we use the following necessary condition for 
LU-equi valence: (pi lj)pij = 17iC/j[(<ti lj)(Tij]C/ ] f 17 i ' 
and therefore tri[(pi lj)pij] = 17jtri[(cri li)au]Uj . 
Since tri[(pi lj)pij] = 1 + a^- XsjUj the equation 



above can be used to determine Vi and Wi (and fc, = 0) 
as those operators which diagonalize tri [(pi li)pu] and 
tri[(cri li)oii] respectively. That is, unless A 3j - = 
Vj Vi and Wi are defined by the equation Vitri[(p\ 
l»)pii]^ f - WJtri[(ai l l )fT ll ]W i t = diag( 7l , 72 ), for 
some 7j. If A 3 ^ = Vi A cannot be orthogonal and 
therefore the unitaries can be determined as explained 
in Sec. IV A. Thus, if p u = U\Ui(JuU\u\ the methods 
described above will lead to the unitaries Vi and Wi in 
Eq. (2) unless pu = p\ 1. 

Another method to compute the unitaries is the fol- 
lowing. Instead of considering the single equation | \I/) = 
Ui . . . U n |$) we use the fact that the basis for the first 
system has been fixed. Therefore, we can use both equa- 
tions given in Eq. (5). Note that the states \^o), l^i) 
and the states \&o), |3>i) are orthogonal respectively. In 
[24] it has been shown that two orthogonal multipartite 
pure states can be perfectly distinguished using local op- 
erations. We are going to use this result now in order 
to determine the unitaries Vi,Wi. It is easy to see that 
for any pair of orthogonal states, |\&o), |^i)i there exist 
local unitaries, Wi such that [24] 



Mi = tr^(|tt ><tfi|) = WWW; 



(6) 



where Ni is a off-diagonal matrix with (Ni)i t 2 = 
ai,{Ni)2,i = bi, for some complex numbers Oi,6j. If 
|<2j| 7^ |6j|, i.e. if NiNj <fc 1, it is easy to see that by im- 
posing the condition that \a%\ > \b%\ this equation deter- 
mines Wi uniquely (up to a phase gate). If |a,| = \b t \ ^ 0, 
Mi is, up to a global phase, a hermitian traceless ma- 
trix. Thus, in this case we would chose Wi such that 
Mj = e la W} DiWi, for some phase a and Di diagonal. 
Defining Vi in the same way for the state |$) we have 
that |\&) ~lu |$) iff Eq. (2) has a solution for the so 
chosen matrices Wi , Vi [34] . 

Before concluding this section we would like to men- 
tion another method to compute the unitaries in Eq. (5) 
for a general state with p\ ^ 1. It is based on the fol- 
lowing observation. Let us denote by |^o) (l^i)) a ( un_ 
normalized) state describing systems 2, . . . , n (2', . . . , n') 
respectively. Then (&~\ u , *o, *i> = iff the state 
I*) = |0) 1*0) + |1) l*i) is either a product state in the 
bipartite splitting system 1 versus the rest, or system i 
versus the rest. This can be easily verified as follows. 
We write |* fc ) - [0), |* fe0 ) + |l),|* fe i), for k e {0,1}. 
Then (tf~ \ u , * , *i) = iff either a) |*oo) |*ii) = and 
|*oi) l*io> = which implies that either |*) = \k) 1 |* fe ) 
or I*) = \k) t (10)! |* fc> + |l)i |*ifc» for k € {0, 1}; or b) 
|*oo) = a l^oi) and |^io) = a\^n), for some a. In this 
case we find |tf) = (a |0) + |1)) 4 (|0) |* i> + |1) |*n»- 
Note that if |^) is a product state in any bipartite split- 
ting, system i versus the rest, then the unitaries Vi,Wi 
and the bit value fcj in Eq. (2) can obviously be easily 
determined. If \^f) is not a product state in this splitting 
then we can combine the two equations in Eq. (5) to 



^(71+72) (g) jj. , 



(*-|..,*o,*i>= (7) 

tf,-,(*-|..,*o,*i>, 



where we used that | \I> ^> = U U \^~) for any uni- 
tary U. This approach will be useful if there are only 
a few unitaries not determined. For instance, if |^) is 
a three-qubit state. Choosing w.l.o.g. i — 2 we have 

(*-| 22 ,* ,*i> = e^+riUsUv (*-\w*o,*i), which 
can then be used to determine C/3, or equivalently V3, W3 
and k^. Of course, the projection onto the singlet state 
can also be performed on more systems. 

In summary, in this subsection we have explained some 
simple ways to compute Vi,Wi and fc, for states which 
have the properties that p, = 1 and that there exists some 
system j such that pj =/= 1. For states with pji ^ pj 1, 
the unitaries can be easily computed using cither that 
(pi T)pu ~lu (ci l)cr 1 j is a necessary condition for 
LU-equivalence or that the states |^o) an d l^i) in Eq. 
(5) are orthogonal. For general states (not requiring that 
Pji =/= pj 1) where there exists a system j with pj =/= 1 
Eq. (7) (and its generalizations) can be used to find new 
conditions on the unitaries. Note that if p\ ^ 1 and 
Pvi = Pi 1 then we have 



I*) 



^|0)(|0)|*oo> + |l)|*oi» 



(8) 



-v/r^p|i)(io> 1*10) + |i) i*ii», 



with (^ij\ ^ki) — l/2SikSji, where p and 1 — p denote the 
eigenvalues of p\. 



V. EXAMPLES 

We are going to employ now the algorithm presented in 
[16] and the results shown in the previous section to char- 
acterize the LU-equivalence classes of up to five qubits. 
We will show that in all these cases it is not necessary 
to determine some unitaries as functions of some others, 
but that it is always possible to determine them directly. 



A. Two- qubit states 

The standard form of a two-qubit state is |W) = 
Ai 1 00) + A2 |11), with Ai > A2 > 0, which coincides with 
the Schmidt decomposition [19]. It is a well-known fact 
that bipartite states are LU-equivalent iff their Schmidt 
coefficients coincide. Let us now demonstrate how this 
result can be rederived with the method presented in [16] 
for two qubits. If Ai 7^ A 2 , i.e. Pi ^ t then, |*) ~ LU |$) 
iff eig(pi) — eig(a\), i.e. iff the Schmidt coefficients Ai 
are the same. For Ai = A2 we have that p\ = pi = 1 and 
therefore the states are LU-equivalent iff eig(p) — eig(a) 
(Lemma 5), which is obviously the case. 



B. Three— qubit states 



following conditions are fulfilled: 



First we transform both states, |\&) and |$) into their 
sorted trace decomposition. If there exists some * such 
that pi 7^ 1, we know that U = Z{on). Without loss 
of generality we assume i = 1. Then we have that the 
states |*) = y/pl\0}\^ ) + ^1 -pi\l) |*i), for some 
P i > 1/ 2 and (ttj I*,) = % and |$) = ^|0) |$o) + 
VI-P2II) |*i), for some p 2 > 1/2 and ($; |$j) = 5 ih 
are LU-equivalent iff 1) eig(pi) = e*<jr(<7i), i.e. iff pi = p 2 , 
and 2) there exist local unitaries U 2 , U 3 and two phases, 
4>, a\ such that 



1*0) = e^t/ 2 ®t/ 3 |$ > 



(9) 



l*i> 



1 U 2 ®U 3 \<$>i 



where 71 = <j>+ai and 72 = (f>—a.\. Since the states are 
not LU-equivalent if p\ 7^ p 2 we assume that p\ = p 2 and 
show now in detail how the unitaries can be computed. 
According to the method summarized in Sec. II we dis- 
tinguish the two cases 1) p\ 2 7^ p\ 1 and 2) pi 2 = p\ <g> 1. 
Since the rank of /)i 2 cannot be larger than two, the sec- 
ond case is only possible if p = 1, i.e. the states |<3/) is 
a product state. Then, the two states are LU-equivalent 
iff the two-qubit states i(0 |\l/) and i(0 |$) are (see Sec 
VA). In the first case we have that either la) at least 
one of the two states |\&j) is not maximally entangled 
or 2b) both are maximally entangled. To investigate the 
case la) we assume without loss of generality that \^o) is 
not maximally entangled and denote by Wi (Vi) the local 
unitaries which map |\&o) (|$o))) hito its standard form 
respectively, i.e. |* > - WM^lOO) + y/T=q[\\l)), 
with gi > 1/2 (|$o) = U 1 V 2 t ( v ^ 2 -|00) + VT — ^ 2 "|ll) with 
q 2 > 1/2). Obviously, |* ) -lu |$o) iff <7i = 0.2- In this 
case, the most general unitaries which transform |$q) 
into |*o), i-e. |*o) = e lll U 2 <g> U 3 |<& ) are of the form 
Ui = WiZ(ai)V i , for some phase a.%. Thus, we have that 
|^I>) ~ LU |$) iff there exists phases Qj such that Eq. (2) is 
fulfilled for ki = fc 2 = k 3 = 0, lq = W-i = 1, and V;, W* as 
defined above. This condition can then be easily checked 
using Lemma 1. In case lb) we have that both, |\&o) 
and l^i) are maximally entangled and are therefore LU- 
equivalent to 1$*). Thus, |\P) is LU-equivalent to |$) in 
this case iff i(fc |4?) is maximally entangled for k = 0,1. 
Since any state in this class is LU-equivalent to the state 
y/p~i\0) |<I> + ) + y/1 — pi |1) l^ - ), the unitaries which map 
two states within this class into each other can be easily 
computed. 

Let us now consider the remaining case where all sin- 
gle qubit reduced states are completely mixed. Since 
pi = P2 = 1 the eigenbasis of pn is maximally entangled 
and therefore LU-equivalent to the Bell-basis (see Sec 
IV A). Thus, any state with pi = 1 Vi is LU-equivalent 
to I*) = |$+) |0) + |$-) |1) = 1 <g> 1 <g> H |*o), where 
|*o) = l/%/2(|000) + |111)) denotes the GHZ-state. 

In summary we obtained the following necessary and 
sufficient condition for LU-equivalence: The two three- 
qubit states |\&) and |$) are LU-equivalent iff one of the 



1) Bi(|*)) = B(|$» = 1 V ? (i.e. Pl = IV*). 

2) There exists some system * such that Ei(\$>)) = 
Ei(\$)) = and Ej(\V)) = Ej(\$)) for some sys- 
tem j =/= i. 

3) There exists some system i such that < Ei(\^>)) < 
1 and Ei(\$>)) = £;(|$)) and either 

3a) Ej(i{k\V)) = Ej(i{k\$}) = 1 for k = 0,1 for 
some system j^i holds; or 

3b) Ejdikl^)) = Ej(i(k\^)) < 1 for one value 
of k € {0, 1} and for the unitaries which can 
be easily and directly determined in this case 
there exists a bit string k and local phase gates 
such that Eq. (2) has a solution. 

For three qubits the polynomial invariants which define 
the different LU-equivalence classes are known [15]. In 
[30] we will compare them to the criterion derived here 
and investigate the measures of entanglement which are 
required to identify the different classes. 

This completes the solution to the LU-equivalence 
problem of three-qubit states. However, in order to il- 
lustrate the method presented in [16], we will apply it 
to the most complicated case, where pi = IV*. We will 
show now that all these states are LU-equivalent with- 
out using the fact that they are LU-equivalent to the 
GHZ-state, |*o)- In other words, we will determine now 
the unitaries {Ui} such that |*) = |* ) = UiU 2 U 3 |<&), 
where |$) = S'jS'^S'g |*). Here, Si are some fixed uni- 
taries. Since the rank of p\ 2 is two, we can compute U 2 
as a function of lq. We find t^tr^G*) (*li |*) {$\)U 2 = 
Ul\i){i\ 2 U 2 = S 2 V /2 (W 1 t |*)<*| 1 Wi|*)('I'|)52, where 
W\ — U\S\. Since the rank of these matrices is one, 
we have (due to the fact that W\ is unitary) that ei- 
ther (0| W\ |0) = (which implies that (1| W\ |1) =0) 
or (lllU/10) = (which implies that (0| W^ |1) = 0). 
Thus, W\ = X kl Z(ai) (up to a global phase), for some 
ki G {0, 1} and some phase a\. Due to the symmetry of 
the state the same holds true for all W, = USj. Thus, 
we have that \^) = UiU 2 U 3 |$) iff there exists ki,k 2 ,k 3 
and phases a, such that |*) = e la ° (g^ Z(on)X ki |W). Of 
course it is straightforward to determine the remaining 
parameters, but in order to continue with the algorithm 
we use Lemma 1 to show for which values of fc, the states 
are up to phase gates local unitary equivalent. The first 
condition, | (i|*)| = | (i| $)| Vi, implies that fci = k 2 = 
k 3 . Then we have \^) is LU-equivalent to |$) iff there 
exist phases a % such that |\I>) = e la ° 0^ Z{a t ) |*), which 
is true iff a = and e »(«i+«2+a 3 ) _ 1 ThuSj us _ 
ing the method above we found Ui = Z(on)X ki Si, with 
e i(a 1 +a 2 +a 3 ) _ j an( j ki = k 2 = k 3 . The unitaries are 
not uniquely defined due to the symmetry of the state. 
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C. Four— qubit states 

In this subsection we will consider the LU-equivalence 
of two four-qubit states, |\&) and |<J>). Similarly to the 
other cases we transform both states, \^f) and |$) into 
their sorted trace decomposition. Like before the solution 
can of course be found using the method summarized in 
Sec. II. However, we will show here that also in this 
case it is possible to determine the unitaries Vi and Wi 
(and the bit values ki) directly. That is, it will not be 
necessary to consider some of the unitaries as variables 
and to determine those unitaries by solving the equations 
which occur in Lemma 1. 

According to the general method, we first distinguish 
the cases 1) there exists some system i with Pi y^ 1 and 
2) pi = 1 for any system i. In the first case we choose 
i = 1 and know that U\ = Z{ot\) for some phase a.\. 
Then we can either have that la) pn ^ p\ 1 or lb) 
Pi2 = pi 1- In case la) W2, V2, and k^ in Eq. (2) can 
easily be determined using the methods presented in Sec. 
IV. Then system 1 and 2 can be measured in the com- 
putational basis leading to four two-qubit states. The 
remaining unitary operators can then be easily found. If 
P\2 — P\ 1 (case lb) we will show next that at least one 
of the unitaries Vi , Wi and the bit-values ki can be deter- 
mined for i £ {3, 4}. First note that the eigenvalues of P34 
arep/2,p/2, (l-p)/2, (l-p)/2 with p ^ 1/2 and there- 
fore P34 is neither 1 nor LU-equivalent to 1— A I* - ) (^ - |, 
for any value of A. Thus, the unitaries Wi, Vi and ki for 
i = 3, 4 can be easily computed unless p 34 =p 3 ®l (or 
P34 = 1 Pa)- However, in this case, since ^34 7^ 1, the 
unitaries V3, W3 and the bit-value k^ (or V4, W4, and the 
bit-value ki resp.) can be easily determined. 

Let us now consider the remaining case where pi = 1 
Vi (case 2). There, all two-qubit reduced states are LU- 
equivalent to 1 + X)i AiSi Si and the reduced states 
are LU-equivalent to each other iff the eigenvalues are 
the same (Lemma 5). We distinguish now the two cases 
2a) P12 7^ 1 and 2b) P12 = 1. In the first case U2 can 
be determined as a function of U\ (see Sec IV). Let us 
apply local unitaries to both states, |<I / ) and |<I>) such that 
Pi2 ; P34 and 012,034 are both Bell-diagonal (see Lemma 
4). We sort the eigenvalues in such a way that if there 
is three-fold degeneracy, then the states are such that 
P12 = P34 = 1-AI*-) (*~| [35]. The resulting states will 
again be denoted by |<3/), |$) respectively. If P12 7^ 1 — 
A \fy~) (^~\, then the unitaries can be easily determined 
(see Sec. IV). In the "worst" case, where p\2 = 1 — 
A|\I/~) (\& _ | we only find U2 — U\. We will show next 
that also in this case U2 , or more precisely V2 , W2 and fe 
in Eq. (2) can be directly computed. That is we will not 
need to compute any of those unitaries as a function of 
some others. 

Using Lemma 4 it is easy to see that any state |\&) 
with P12 = 1 — A |\I> _ ) (it - 1, is LU-cquivalcnt (up to 
a global phase) to a sta te |$+) |$+) + e 471 |$") |$") + 
|*+) |*+) + Vl - Ae* 73 |* _ ) |* _ ), for some phases 



change the sign of two states out of the four Bell states, 
we can choose 71, 72 < 7r. We are going to show next that 
two states of this form with the choice 71 , 72 < tt are LU- 
equivalent iff the complex coefficients which occur here 
coincide. 

Let us denote by U m b the 4x4 unitary matrix, which 
transforms the computational basis into the magic ba- 
sis, i.e. U mb |00) = |$+) , U mh |01) = -i |$-) , U mb 1 10) = 
I* - ) , U mb |H) = —i |* + )- It is a well-known fact that 
for any U\, U2 unitary we have that U mb U\®U2U m b = O, 
where O is a real and orthogonal 4x4 matrix. Further- 
more, it is easy to see that Oi = U mb Ui UiU m b can be 
written as O l = O l ® |01) (01| (U mb |01) = |^ - )), where 
Oi is a three dimensional rotation. 

Since p\2 = v\2 = 1 — A \^~) (^~ | and therefore P34 = 
(T34 = 1 A|\&~) (fy~\ we know that there exist local 
unitaries, Ui which map |$) into |\I/) iff U2 = U\ and 
Ui = U3. Thus, we have |\fr) ~lc/ W) iff there exist real 
orthogonal matrices, Oi — Oi © |01) (01 1, with Oi is a 
three-dimensional rotation such that 



^') = ul b ®ul b 



m = Oi O3 



$ 



(10) 



Note that 

e* 72 111) 111) + 






= 1 00) 1 00) - e* 71 1 10) 1 10) - 
"Ae* 73 |01) |01), and that the phases 



$) = u mb < 



7j are not local phases. Similarly we have 

U mb |$) = e* 7o (|00) |00) + e' 71 |10) |10) + e* 72 |11) |11) + 

V 1 — Ae* 73 |01) |01)), for some phases 7, and coefficient 
A. Using now that the real and orthogonal matrices 
Oi,Oz are of the form Oi © |01)(01| it is easy to see 
that I*) ~ LU |$), i.e. Eq. (10 ) is sa tisfied iff {e i7 *}f =1 = 

{e i7 *}f = i and VT^Xe^ 3 = y/l-Xe^' 3 . 

For the case 2b), where p\2 = 1, which implies that 
p 34 = 1 we write |*) = £V \hj) \ipi,j) witri d'ijli'kl) = 
dikdji- Since these states form an ON-basis we can 
find a 4 x 4 unitary U such that \ipij) = U\ij). Re- 
call that any two-qubit unitary operator U can be writ- 
ten as U = U\ U2UdV\ V2, where Ud, the non 
local content of U, is diagonal in the magic basis, i.e. 



U d 



a i(<p 1 X®X+4> 2 Y®Y+4> 3 Zt»Z) 



for some phases 



o«72 



[23]. Note that Ud can be made unique by imposing 
certain conditions on the phases, 4n [25]. We trans- 
form the state by local unitary operations into the form 
I12 U d Y,ij \ij) \%j) = I12 U d |$+) 13 |$+) 24 . Then 
the two states, |*) = 1 12 U d ($) ^ l$ + )i 3 l* + ) 2 4 and 
|$) = I12 U d ($)J2ij I* + >13 l $+ ) 24 are LU-equivalent 
iff Ud(&) = Ud(&)- Thus, these LU-equivalence classes 
are characterized by Eij(\^>)) — 2 for some systems i,j 
and the three parameters, 4>\,4>2,4>3i which define the 
non-local content of Ud- In Sec VIII we will give a phys- 
ical meaning to these parameters and discuss its gener- 
alization. 

Using the fact that any four-qubit state with p\2 = 1 
is LU-equivalent to the state |\P) = li2 0^d^i 7 - l*j) l*J')> 



7i. Since the operations Sj Ej for « e {1,2,3} always where C/ d = t/ m &diag(l, e 2 * 11 , e 



i0l P i02 o*<feurt 



^mb' for SOme 
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phases 4>i , it is also easy to rederive the result that there 
exists no four qubit state with the property that all two- 
qubit reduced states are completely mixed, i.e. p^ = 1 
Vi,j- This can be seen as follows. Since |\I/) = li 2 
UdEu ly> l*i>. wefindp 13 = tr 4 (t/ d |$+) (*+ 



i:r 



)tiU\) 



Then the conditions p\$ = p23 = 1 imply that cos(0j) = 



and cos( 



= Vi, j. Since it is impossible to fulfill 



those equations simultaneously, we have that there exists 
no four-qubit state such that p^ = 1 Vi, j. This implies 
that the case 2b) is actually contained in 2a). In fact, it 
corresponds to the case where A = 0. 

In summary, for the four qubit case we have the fol- 
lowing possibilities: 

la ) There exist some systems i and j with Pi ^ 1 and 
Pij 7^ Pi 1: Without loss of generality we choose 
% = 1 and j = 2. Then Wi = Vi = 1 and fci = in 
Eq. (2) and V2, W2 and fe can be easily computed 
using the methods presented in Sec. IV. 

lb) There exists some system i with pi ^ 1 and some 
system j with p^j = pi ® 1: Without loss of gen- 
erality we chose i = 1 and j = 2. Then the uni- 
taries Wj, Vi and fcj for i = 3, 4 can be easily com- 
puted by considering ^34 which can neither be 1 
nor 1 — A \*&~) ( i &~\, for any value of A. 

In both cases at least for two systems the operators 
Vi,Wi and the bit values fc, can be determined. Thus, 
measuring those two systems in the computational basis 
leads to four equations for two-qubit states. The missing 
operators Vi , Wi and the bit values ki can then be easily 
computed. The states are LU-equivalcnt iff there exist 
some phases c^ such that Eq. (2) has a solution, which 
can easily be checked using Lemma 1. 

2a) For any system i, pi = 1 and there exists a sys- 
tem j such that pij =/= 1, for some system i: With- 
out loss of generality we chose i = 1 and j = 2. 
First we apply local unitaries to both states, |\&) 
and |$) such that pi2,P34 and (712,034 are all 
Bell-diagonal (see Sec IV A). The resulting states 
will again be denoted by |\P), |$) respectively. If 
P12 7^ 1 - A |*~) (\ff~| for some A e R (which im- 
plies that also that ,034 7^ 1 — A |$~) (*&~\ for some 
A e R), the unitaries Vj,Wj and fc, in Eq. (2) 
can be directly computed using the methods of Sec. 
IV. If P12 = 1 - A |*~) (*~| for some A e R we 
map I*) into the form |$+) |$+)+e 471 |$~) |$~) + 
e n2 |^r+^ |$+)+VAe i T 3 I*") I*"), for some phases 
7i with 71.2 < 7r. Two states of this form are LU- 
equivalent iff their complex coefficients which occur 
here coincide. 

2b) For any system i, p h = 1 and there exists a sys- 
tem j such that p^ — 1, for some system i: 
Without loss of generality we chose i = 1 and 
j = 2. In this case the state is LU-equivalent 
to the state |*) = 1 12 (8) ftdWX^- Ki) \ij)> where 



£/<*(*) = C/ mh diag(l,e*,e^,e^)C/^ b canbecho- 
sen uniquely. Then two states are LU-equivalent 
iff Ud(&) = Ud(<&). Note that since there exists 
no four-qubit state with pij = 1 Vi, j case 2b) is 
contained in 2a). 

Similarly to the three-qubit case, we will consider 
now the most complicated case (for the algorithm pro- 
posed in [16]) and show how the unitaries which trans- 
form two LU-equivalent states into each other can be 
determined. The most complicated case for four qubits 
is the one where p^ = 1, for some systems i,j which 
we chose to be 1,2. Thus, we consider the example 
I*) = I12 C/34 |$ + ) 13 | < E >+ ) 2 4, where we choose t/34 = 
L/ m6 diag(l, e^, e^, l)L/^ 6 such that p 12 = p 3 4 = P23 = 1- 
It can be easily shown that pi% = p 2 4 = 1/4(1 + X <£> 
X + cos(4>)(Z Z — Y Y). Our aim is to determine Ui 
such that I*) = Ux . . . Ui |$), where |$) = S\ . . . S 1 ] |#), 
for some given unitaries Si. Since P13 7^ 1 and P24 =/= 1 
we can compute U3 (U4) as a function of U\ (U2) re- 
spectively Considering Eqs. (4) for all values of I and m 
simultaneously we have P13 = l/4(l+AT®AT+cos(0)(Z<8> 
Z -Y F)) = C/iC/ 3 (5' 1 t S' 1 t pi35'iS , 3)t/ 1 t t/3 t . It is straight 
forward to see that the last equation can only be fulfilled 
if U X S\ = U 3 Sl = E fe , for k <= {0, 1,2,3}, where S = 1. 
Similarly we find U2SI — UiS\ — £;. Thus, we have 
!7j = SfcjiSj, where k\ ~ k 3 = k and /c 2 = &4 = ^- It is 
straightforward to show that |\I/) = SfeS;SfcE/ |^) forcer- 
tain values of fc, Z (e.g. fc = 0, 1 = 3 or fc = 1, 1 € {0, 1}). 
Again, the reason, why the unitaries are not uniquely de- 
termined using this method is because of the symmetries 
of the state. 



D. Five— qubit states 

Instead of considering now, similarly to the other cases, 
all possible classes of five-qubit states we consider here 
one of the hardest examples to illustrate the method pre- 
sented in [16]. First, we are going to construct a five- 
qubit state, |\&), which has the property that all the 
two-qubit reduced states are completely mixed. Then, 
we consider the two states I*) and |i>) = Si® . . .®Ss \^f) 
for some local unitaries Si and compute the unitaries, Ui, 
which map |$) into |\&), using the algorithm presented in 
[16] and summarized in Sec III. We will show that also 
in this case it will not be necessary to determine any of 
the unitaries as a function of some others, but that it will 
be possible to determine them directly. 

In order to construct the 5-qubit state, |\&) with 
Pij = 1 Vi,j we write |tf) = |0) |* > + |1> |*i>, 
where (*, |\I/j) = 5ij. As shown above, any four- 
qubit state which has the property that P12 = 1 is 
LU-equivalent to a state I12 UdJ2i \ij) \ij), where 



Ud = t/ m fcdiag(l,e Ml ,e M2 ,e M3 )£/ r T n6 . Imposing now also 
that P23 = P14 = 1 and that the phases on fulfill 
< at < it, we find |*o) = I12 U dl J2 itj \ij)\ij), 
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where U dl = U mb diag(l,e ia \e ia \ 1)17^ 



It is easy to 



li2 0f/d 2 Ei, j Ni)Ni).with 
M1 1 ) ?7 m6- xt is straight for- 



see that two states, |\&o) , |*i), of this form are orthogo- 
nal to each other iff a 2 = a\ + it. We consider now the 
5 qubit state [*) = |0) |* ) + |1) |*i) with |* ) = I12 <S 
Udi Eij l*J> l*i> and 1*1 
t/d 2 = t7 m6 diag(l,-e iai ,-e 
ward to show that all two-qubit states, p^ are completely 
mixed. Note that |*) = |+)(|$+,$+) + | *+,*+)) + 

e i«|_)(|$- ) $-) + |*-,*-)). 

Let us now consider the state | $) = 5i . . . ® S$ | *) 
for some local unitaries Si and compute the uni- 
taries, Ui, which map |4>) into I*). Since pij = 1 
we choose U\,U 2 as parameters. Note that p\ 23 7^ 1, 
since it can have at most rank 4. Thus, we can 
compute the unitary U 3 as a function of U\ and C/2 
considering P123, which can be easily shown to be 
1/8(1 + X ® X ® X). We have [7* |#) = U 1 U 2 lU i U 5 |$) 
and therefore E^tr^d^)^ <*il I*) (^1)^3 

tr^r/M | W)i2 (ij | ^^ | $) ($!) for any ^ fcj L It 

can be easily shown that tr-, 3 (|00) 12 (11| |*) (*|) = X/2 

and that tr-^M I 00 )i2 < n l U \U 2 I*) ($1) = xS 3 XSl 
where x depends on U\,U 2 . Since only x depends on U\ 
and U 2 , U 3 can be directly computed (not only as a func- 
tion of U u U 2 ). We find U 3 - ° ia * x < ' 



Thus, denoting by 



' ■"*" S^ for some phase a 3 . 



* 



Hz I*) we have I*) 



-LU 



m 



iff there exist local unitaries U\,U 2 ,U&, U5 and a phase 



a 3 such that 



*\ = U 1 U 2 Z(a 3 )HSlU 4 U 5 |$), where 
53 is determined. Projecting now the third system 
onto |0) we find a state of system 1245 with the 
property that p 2 4 oc 1 + X X . Imposing then 
the necessary condition p 2i = U 2 JJ^Q 2 JU\lj\ of LU- 



equivalence leads immediately to U 2 



Ja 2 X ct 



S\ and 



U4, = e laiX S\. Similarly we find the other unitaries. 
Thus, we have |*) ~lu |$) iff there exist phases a,i such 

that 0f = i H i l*> = eM ° <S>. ^(«i)(®?=i H i S l) l $ >- Thc 
existence of the phases can be easily verify either by 

looking at the coefficients in the computational basis, or 

by employing Lemma 1. 



VI. LOCC INCOMPARABILITY 

The results presented here lead also to conditions for 
the existence of more general operations transforming one 
state into the other, namely Local Operations and Clas- 
sical Communication (LOCC). This is due to the fact 
that two multipartite states, having the same marginal 
one-party entropies, are either LU-equivalent, or LOCC- 
incomparable [26, 27], i.e. none of the states can be 
mapped into the other by LOCC. 

In this section we will show that for any n > 2 there 
exists a pair of n-qubit states, {|\&) , |<I>)} such that for 
any bipartition A/B, where A contains a qubits and B 
n — a qubits, eig(pA) — eig(&A), but the states are not 
LU-equivalcnt. In particular, the entropies of the re- 



duced states of any subsystem coincide, i.e. all bipartite 
entanglement, measured with the von Neumann entropy 
of the reduced states is the same for both states. Since 
the eigenvalues of all single qubit reduced states coin- 
cide those states are not even LOCC comparable [26, 27]. 
Surprisingly, in those examples |#) will be the complex 
conjugation (in the computational basis) of |*). That is, 
for any n (n > 2) there exist n-qubit states which are 
not even LOCC comparable to its complex conjugate. 
The consequence of the existence of these states will be 
discussed in Sec VIII. 

Note that, for 3-qubit states examples of such states 
have already been presented in [28]. There, the fact that 
the states are not LU-equivalent has been proven by em- 
ploying a polynomial invariant of degree 12. Here, we 
will use the necessary and sufficient condition for LU- 
cquivalence presented in [16] and Sec. Ill to prove that 
thc considered n qubit states are not LU-equivalent. 
First we will present a three-qubit state, |\fr) which is 
not LU-equivalent and therefore not even LOCC compa- 
rable to its complex conjugate. Then we will generalize 
this example to n-qubit states. 

We consider the LME states |W) = U 123 (<t>) \+) m and 
1$) = U 123 (0+7T ) \+)® 3 , where t7i 23 (a) is the three-qubit 
phase gate defined by U 123 = 1 - (1 - e 1 *) |111) (111|. 
As mentioned in Sec III an arbitrary LMES, |\E f ), can be 
easily transformed into its trace decomposition, |*t r ), by 
applying the local unitary operations HZ(cf>i), where <pi is 
chosen such that cot(<fii) — Tpr- For the symmetric state 
I*) we find (Xi) = 1/4(3 + cos((f>)) and (Fj) = sin(0)/4 
and therefore cot(4>i) = cot(</>) + 3csc(</>) for i = 1,2,3. 
For <j> — 7r/2 the marginal entropies of I*) and [<&), which 
is equivalent to the complex conjugate of I*) in this case, 
coincide. However, it is easy to show that \^tr) I ■ \^tr)* 
is not a product state and therefore the states I*) and 
I*)* are not LU-equivalent (see Lemma 3). Moreover, 
due to the fact that the eigenvalues of all the reduced 
states are the same for [*) and |**), those two states are 
not even LOCC comparable. Note that those two states 
have the same bipartite entanglement (considering any 
bipartite splitting) and the same value for the tangle [5] , 
the value of which is the same for a state and its complex 
conjugate. 

Let us now generalize this example to n-qubit 
states (for n > 2). That is the two n-qubit states, 
I*) = E7i,...,„(7r/2) |+) " and |$) = |**) have the 
property that eig(pA) — eig{<?A) for any subsystem A. 
However, the states are not even LOCC comparable. In 
order to prove that we first note that the eigenvalues of 
Pa and p A coincide for any subsystem A. Furthermore, 
since the state is symmetric with respect to particle 
exchange, all single qubit reduced states coincide. 
They are of thc form p = |+) (+| + 2-"( v / 2((e ra - 
1) |1) (+| + {e~ ia - 1) |+) (1|) + \e ia - 1 | 2 |1) (1|), with 
eigenvalues 1/2(1 ± 2~V8 - 2 2 +™ + 2 2 ™). Thus, none 
of the reduced states is proportional to the identity 
and therefore the states are LU-equivalent iff their 
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standard forms coincide. Since those states are LMESs, 
we know that their trace decompositions are of the form 
[HZ(a)]® n |\&), where a is determined via the equation 
cot(a) = -Wfj = 1 — 2™ _1 . It is straightforward to show 

that \% r ) = e M "/ 2 [cos(a/2) |0) - isin(a/2) |l)]® n + (i- 
l)e m ™/ 2 [|0) — |l)]® n and therefore none of the coefficients 
in the computational basis vanishes. There exist now 
several ways to prove that |^) is not LU-equivalent to 
its complex conjugate. We could either compute the 
standard form, \^f s ), and show that it does not coincide 
with the one of |^*), i.e. show that |\& s ) is not real, or 
we could employ Lemma 1 or Lemma 3. We will use here 
Lemma 3 to show that the sorted trace decompositions 
are not related to each other by local phase gates, which 
proves that the states are not LU-equivalent. Since 
the hrst condition (condition (i)), | (i| * tr )| = | (i| $* r )| 
is obviously fulfilled, we have that |*) ~ LU |v]>*) iff 
|*) /. |**) is a product state. Since |\P) is symmetric, 
this last condition is fulfilled iff there exists a single qubit 
state \<f>) such that |*) /. |**) = |<£)® n . In other words, 
the states are LU-equivalent iff there exists two phases ao 
and ai such that |* tr ) = e la °Z( ai )® n |** r ). This last 
equation is fulfilled iff there exists a phase a\ such that 
tfi,..,n(7r/2) \+f n = e ia °V® n U h ...,„(-7r/2) \+f n , 
where V ai = Z(-a)HZ(ai)HZ(-a) = 

cos(ai)Z(— 2a) + is'm(ai)e~ la X. Rewriting this 

condition we have that the states are LU-equivalent iff 
there exists a phase a\ such that |+) n + (i — 1) |1) " = 
e iao ((V ai |+»® n + (-i - l)(V ai |1))®"). It can be 
easily shown that this condition can only be ful- 
filled if V ai |+) ex |ei) and V ai |+) ex |e 2 ), where 
\ e i) G {|1) ) 1+)}) f° r i = 1)2 and ei c^ e 2 . We consider 
the two possible cases, (a) V ai |+) = a\l) for some 
a e QJ and (b) V ai |+) = a\+) for some a G QJ. Case 
(a) is possible iff cos(ai) + is'm(ai)e~ la — and 
cos(ai)e~ 2M + isin(ai)e~ M = a. It is easy to see 
that the first condition cannot be fulfilled since a is 
determined as mentioned above. Since V ai |+) = a\+) 
(case (b)) implies that cos(ai) = it is also easy to 
see that in this case V\l) = ±ie~ ta |0) and therefore 
V|l) ^ a\a) with \e t ) G {|1),|+)}, for i = 1,2. This 
proves that for any n > 2, the two n-qubit states, |*) 
and \ty*) are not LU-equivalent. 



VII. LU-EQUIVALENCE OF MIXED STATES 
AND D-LEVEL SYSTEMS 

We will show here that the criterion of LU-equivalence 
presented in [16] serves also as a criterion of LU- 
equivalence for certain mixed and also for certain mul- 
tipartite states which describe a system composed out of 
d-level systems. 

For instance, if we want to find out whether or not two 
mixed states are related to each other by local unitarics 
and if there exists at least one eigenvalue of the mixed 
state which is not degenerate then the same method can 
be used. This is due to the fact that the unitarics can- 



not change the eigenvalues, thus, if p ~lu Q then it must 
hold that |*) ~ LU |$), where |*) (|3>)) denote the eigen- 
states to the non-degenerate eigenvalue of p, (a) respec- 
tively. In order to check then if the two mixed states are 
LU-equivalent, one first uses the algorithm to determine 
the local unitarics which transform |3>) into |\I/). Those 
unitaries must also transform a into p, which can then 
be easily checked. 

The criterion for LU equivalence for pure states can 
also be employed for mixed states if there does not exist 
a non-degenerate eigenvalue, but one which is two-fold 
degenerate. Let us denote by |\&o) > l^i) an d |3>o) > |$i)j 
the eigenvectors corresponding to the two-fold degener- 
ate eigenvalue of p, a respectively. As before we have 
that if p ~lu C) then there exist local unitaries, Ui 
such that |*o)(*o| + |*i)(*i| - ®i C7i(|* > <*o| + 
|$i) ($i|)<8>i UJ. This equation is fulfilled iff there ex- 
ists a 2 x 2 unitary, V, with |* fe ) = ® 4 U t Y,i Vkl \®i) , for 
k = 1,2. Note that this is equivalent to finding the local 
unitaries which map the state |$) = |0) |$ ) + |1) |$i) 
into |*) = |0) |* ) + |1> |*i)- Thus, solving the LU- 
equivalence problem of mixed n-qubit states, where one 
eigenvalue is two-fold degenerate is equivalent to solving 
the LU-equivalence problem of n + 1 qubit states, where 
Pi = 1. ' 

Suppose now that p is a n-qubit mixed state and its 
eigenvalue with the smallest degeneracy is ^-fold degen- 
erate. We denote by \$>k) (l^fe)) the eigenstates of p, 
(a) correspond to this eigenvalue. Then p c^lu & implies 
that J2k\ k )\^k) = V<g> t U t J2 k \k)\<P k ), where V is a 
k x k unitary and all the other Ui are single qubit uni- 
tary operations. The idea is then to first fix the unitarics 
Ui using the algorithm presented in [16] and at the end 
try to fix the unitary V. 

The LU-equivalence problem for d-level systems can 
be investigated in a similar way. However, due to the 
additional degeneracy which can occur in this case, the 
situation gets more complicated. For instance, if a state 
|\&) describes a system composed out of <i-level systems, 
then pi, which describes a single <i-level system can be 
I -fold degenerate, where I < d. Thus, in this case the 
unitaries occurring in Eq. (2) can in general not be de- 
termined up to local phase gates. If there is no degener- 
acy similar methods can of course be applied to solve the 
problem of LU-cquivalence. 



VIII. MULTIPARTITE ENTANGLEMENT 

The algorithm presented in [16] cannot only be used 
to solve the LU-equivalence problem, but allows us also 
to gain a new insight into the entanglement properties 
of multipartite states. Within the algorithm the classes 



Pn 1 ,...ni,k 



¥= Pn u ...n, 1 and p ni . 



■ ni.k 



for some subset of qubits, rii, . . .ni,k, are distinguished. 
In the first case the unitary Wk can be computed as a 
function of the unitaries U ni , . . . U ni , whereas it cannot 
(using the proposed algorithm) be computed in the sec- 
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ond case. Therefore, in this case a new variable, Uk is re- 
quired. As explained in [16] those classes correspond also 
to different entanglement classes. For instance, applying 
any von Neumann measurement on the first subsystem 
described by the state |\1/), with p\2 = pi 1, always re- 
sults in a state where the second system is maximally en- 
tangled with the remaining systems, independent of the 
measurement outcome. Obviously this is not the case 
for a state with pn ^ p\ 1. This suggests that, in 
order to understand how a many-body system can be 
entangled, one first identifies the entanglement class (as 
described above) to which the state belongs to. Note 
that this classification is based on multipartite, not bi- 
partite entanglement properties. However, it is easy to 
perform this classification since one only needs to con- 
sider the reduced state of certain subsystems. Within 
the identified entanglement class it is then feasible to un- 
derstand how multipartite entanglement can be qualified 
and even quantified. For instance, as we have seen in 
Sec. V, the LU-equivalence classes of four-qubit states 
with pij = 1, for some systems i and j are characterized 
by three parameters. This is due to the fact that any 
state in this class (choosing w. 1. o. g. i = 1, j = 2) 
can be written as |*) = 1 12 Ud |<J? + ) 13 |$ + ) 24 , where 

U d = e «.<l>iX®X+<t> 2 Y®Y+<t> 3 Z®Z) ; for some phases fa 

Thus, also the entanglement contained in such a state 
is completely characterized by E^d^)) = 2 and the 
three phases, 4>i. Recall that any two-qubit gate, U, 
can be decomposed as U = U\ UiUdYx V2, with 
Ud as above denotes the non-local content of the gate 
U. Using all that allows us to give the three parameters 
4>i the following physical meaning. Recall that the state 
I*) = I120C/34 |$ + ) |$ + ) is the Choi-Jamiolkowski state 
corresponding to the operation U [17]. That is, given 
the state |\I/) the operation U can be implemented using 
just local operations [36]. This shows that the non-local 
properties of a four-qubit state for which there exists a 
maximally entangled bipartite splitting between two ver- 
sus two qubits is completely characterized by the amount 
of entanglement which can be generated using this state 
as the only non-local resource. 

Of course, this new insight in characterizing multipar- 
tite entanglement by the amount of entanglement which 
can be generated using this state as the only non-local 
resource can be generalized to an arbitrary state, inde- 
pendent of the dimension and even for mixed states. Note 
that the quantum operation corresponding to a state de- 
scribing n subsystems is acting only on \n/2~\ systems. 
Like in the example of four-qubit states, the correspond- 
ing operation is acting on two qubits. This fact simplifies 
the characterization of multipartite entanglement, since 
e.g. the non-local properties of two-qubit operations are 
very well understood. It should be further noted here 
that the operation corresponding to the state, |\&) via the 
Choi-Jamiolkowski isomorphism is unitary iff the state 
has the property that it is maximally entangled in the 
considered bipartite splitting. 

Let us point out here that the algorithm gets more and 



more complicated the larger the number of systems, I is 
for which p ni .... ni = 1 for any choice of Hi, . . . ni, since 
then I unitaries have to be considered as variables. In 
the worst case, where any bipartition of |~n/2] qubits is 
maximally entangled with the rest, |~n/2] unitaries have 
to be considered as variables. It is known however, that 
only for very few values of n such states exist [29] . On the 
other hand, the more systems are maximally entangled 
with the rest, the less parameters remain to characterize 
the LU-equivalence class. Like, for instance in the exam- 
ple of four qubits, the class with p^j = 1, for some systems 
i,j can be characterized with only three parameters. For 
the other extreme case of generic states all the parame- 
ters occurring in the standard form determine, like in the 
bipartite case, the entanglement contained in the state. 

Another important insight into multipartite entangle- 
ment which we gained here is the fact that for any n 
there exists a n-qubit state, \^f), which is not LOCC 
comparable to its complex conjugate (Sec VI). Thus, 
the non-local properties of |\P) and |\£*) seem to be 
really different. Since the mapping |W) — > |\I/*) corre- 
sponds to the redefinition of the complex unit i by —i, 
one might expect that this change does not lead to any 
new physics. In fact, for any observable O, we have 
(*| O I*) = (**| O* I**). Thus, whatever measurement 
outcome we can get by measuring a system described 
by the state |W), the same outcome can be obtained by 
measuring O* on |^*). Due to that, there will not ex- 
ist a physical measure which is capable of distinguishing 
those two states. This shows that it will not be possible 
to characterize all LU-equivalence classes by entangle- 
ment measures. 

This suggest the introduction of a function, I\, with 
Ji(|*)) = if I*) ~ LX j I**) and 7i(|*» = 1 otherwise. 
If 7i(|\&)) = 1 the Hilbcrt space should be divided into 
two subsets, one containing |\&) and the other containing 
l^*). After making this distinction one proceeds inves- 
tigating the non-local properties of the state \^>) within 
the subset associated to it. 

In [30] we will follow the approach to investigate the 
multipartite entanglement properties in the way outlined 
here. In particular, we will consider multipartite state, 
describing several qubits, and will introduce the func- 
tion which determines if a state is LU-equivalcnt to its 
complex conjugate or not. Moreover, we will analyze the 
entanglement contained in the state by investigating the 
amount of entanglement which can be generated using 
this state as the only non-local resource to implement 
quantum operations on a smaller system. 



IX. CONCLUSION 

We used the criterion of LU-equivalence of multipar- 
tite pure states to derive the different LU-equivalence 
classes of up to four qubits. For five-qubit states, which 
can be treated analogously, it is shown that the most 
complicated class of states, where all two-qubit reduced 
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states are completely mixed can be easily considered us- 
ing the algorithm developed in [16]. Even though it is in 
principle necessary to determine some of the local uni- 
tary operations as a function of some others, it is shown 
that for those cases this is, in fact, not required. That 
is, the unitaries can always be directly computed. The 
algorithm suggests to distinguish different classes of en- 
tangled states, like the one where p\2 ^ p\ ® 1 and where 
P12 — Pi ® 1- We considered here all the possible classes 
and showed that within certain classes new, operational 
entanglement parameters can be identified which com- 
pletely characterize the non-local properties of the states. 
For instance, it has been shown that any four-qubit 
state for which one two-qubit reduced state is completely 
mixed is LU-equivalent to a state |*) = \®Ud |<& + ) [$+), 
where U d = e ^iX®x+ a2 Y®Y+ a3 z®z) with a . £ R ig the 

non-local content of a two-qubit gate [23]. The state |\&) 
is the Choi-Jamiolkowski state corresponding to the op- 
eration Ud [17]. Thus, Ud can be implemented by local 
operations if the state |\&) is used as a resource. This 
new approach of characterizing the entanglement of a 
multipartite state by the entangling capability of the op- 
eration which can be implemented using the state as the 
only non-local resource can be generalized to arbitrary 



states. Moreover, we derived examples of n-qubit states 
(for n > 2) which arc not LOCC comparable to their 
complex conjugate. This observation suggests the intro- 
duction of a new measure, which distinguishes the cases 
|*) ~lu |**> and |*) ^ LU |**). If the states arc not 
LU-equivalent, two different subsets of the Hilbert space 
should be considered, one for |\P) and one for |^*), in or- 
der to further investigate the properties of multipartite 
entangled states. In [30] we will prove that the examples 
of n-qubit states, |\P), which are presented here, cannot 
even be mapped into their complex conjugate by allowing 
stochastic LOCC (SLOCC). That is, it is not possible to 
transform \^>) into |\I/*) by local operations even with an 
arbitrary small probability of success. 
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